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The cosmic ray interaction event generator Sibyll is widely used in extensive air shower sim- 
ulations for cosmic ray and neutrino experiments. Charmed particle production has been added 
to the Monte Carlo with a phenomenological, non-perturbative model that properly accounts for 
charm production in the forward direction. As prompt decays of charm can become a significant 
background for neutrino detection, proper simulation of charmed particles is very important. We 
compare charmed meson and baryon production to accelerator data. 



I. INTRODUCTION 

The cosmic ray event generator Sibyll is a rel- 
atively simple model that is able to reproduce many 
features of hadronic interactions in fixed target and 
collider experiments. Sibyll is based on the dual par- 
ton model (DPM) [i-lH, the Lund Monte Carlo algo- 
rithms [1, and the minijet model [8l-[ll|. allowing 
multiple hard and soft interactions. The hard and soft 
interactions are discriminated by an energy-dependent 
transverse momentum cutoff. The hard interaction 
cross section is calculated according to the minijet 
model. Diffraction dissociation is implemented as a 
two-channel eikonal model based on the Good- Walker 
model [T2j . For hadron- nucleus interactions, the inter- 
action probability for each nucleon inside the nucleus 
is calculated based on the impact parameter distri- 
bution. The total interaction cross section is calcu- 
lated using the Glauber scattering theory [Ts'l. For 
a nucleus-nucleus interaction the semi-superposition 
model [lil is used to determine the point of first in- 
teraction for the nucleons of the projectile nucleus. 
The fragmentation region is emphasized as appropri- 
ate for air shower simulations. Version 2.1 has been 
available since 1999. 

The new version 2.2c has charmed quark and anti- 
quark added to the u, d, s quarks and anti-quarks, 
and gluons which are present in version 2.1. Small 
updates have also been made to give better agree- 
ment with experimental data. Addition of charmed 
quark is especially relevant to air shower simulation 
studies of neutrino detection: the short-lived charmed 
particles create "prompt" muons and neutrinos in the 
atmosphere with a harder spectrum and isotropic dis- 
tribution, whose contributions are searched for in the 
atmospheric muon and neutrino spectrum. A proper 
representation of charmed particles is thus very impor- 
tant for simulations involving neutrino experiments. 

In the following sections, we describe how the 
charmed quark is added to Sibyll, and compare the 
charmed meson and baryon predictions with exper- 
imental data. Hereafter, a "g quark" includes both 



the quark and anti-quark unless specified. 



II. CHARM INCLUSION 

The following charmed particles have been added: 
mesons D, D* families, rjc, and J/ip] and baryons 
Ec, S;, Sc, S* families, A+, and fl°. 

Charmed particles are formed throughout the whole 
fragmentation process. The basic assumption is that c 
quarks form in all processes where s quarks form, but 
with lower frequency. This is realised in the Monte 
Carlo by replacing production of an ss pair by a cc 
pair with a frequency Pc/Ps = 0.004. The numerical 
value is adjusted for a best fit to the data. For com- 
parison, Ps/Pd = Pus /Pud = 0.3. Valence s-quarks 
(for example in kaon projectiles) are excluded from 
this substitution process. This procedure ensures that 
all charm production coming from hard and soft pro- 
cesses are covered. The cross section is controlled by 
the branching ratio. 

In a collision, the leading particle carries a valance 
quark or diquark of the projectile, in contrast to the 
non-leading ones. Leading particles therefore are more 
abundantly produced in the forward region, i.e. large 
Feynman xp, resulting in an asymmetry between lead- 
ing and non-leading particles. This leading quark ef- 
fect is well established in charm hadroproduction. For 
example, the quark composition of a proton projectile 
favours A+ over A^ as leading particles. Charmed 
particles formed via fragmentation allows the leading 
quark effect to be observed 15] . 

The momenta assigned to the charmed particles dif- 
fers from non-charmed particles to give a harder spec- 
trum. While the non-charmed particles follow the 
Lund fragmentation function, the heavier charmed 
particle fragments according to the Peterson/SLAC 
fragmentation function jT^ 
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where z is the fraction of the new particle energy with 
respect to the parent quark or diquark, and the free 
parameter eg, which is inversely proportional to the 
square of the heavy quark mass, has been set to 2. The 
primordial pt assigned to the newly formed c quark 
or diquark pairs during fragmentation has a Gaussian 
distribution with an energy-dependent mean value 



10' 



Po 



0.08 logi 



30 GeV 



GeV/c, (2) 



where pq = 1.5 GeV/c and 1.0 GeV/c for baryons and 
mesons, respectively. This setting results in a higher 
(pt) for charmed particles, compared to other quarks 
and diquarks where the (pt) ranges from 0.3 to 0.6 
GeV/c. 

Charmed particles are produced in the interactions 
of high-energy primary cosmic rays. The quick decay 
of these charmed particles form the "prompt" muons 
and neutrinos in the atmosphere. These prompt lep- 
tons are expected to have an isotropic angular distri- 
bution and a harder spectrum than particles produced 
from ordinary decay of pions and kaons. Because of 
their relatively hard spectrum and isotropy, prompt 
neutrinos are a significant background in the search 
for astrophysical neutrino signals above 100 TeV. 

Additional minor modifications with respect to 
SiBYLL 2.1 have been made in version 2.2c. These 
include improvements in treatment of diffraction, in- 
creasing the s quark fraction and fixing a bug in the 
energy dependence of the cutoff of pt- These changes 
give a higher multiplicity and a smoother transition 
to diffraction in the forward direction. 



III. COMPARISON WITH DATA 
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FIG. 1: Differential cross section distributions of xf (up- 
per panel) and (lower panel) for p-p collisions producing 
Ds from Sibyll 2.2c plotted against data from LEBC- 
EHS {Eiat = 400 GeV) and LEBC-MPS {Eiat = 800 
GeV). For clarity purpose, the simulation results and data 
for LEBC-MPS has been multiplied by factor 10. 



Most information on differential behaviour is ob- 
tained from fixed target experiments, with hadron- 
hadron or hadron-nuclei collisions. Distributions of 
Feynman xp and transverse momenta p^ from the 
LEBC bubble chambers, E769, and SELEX experi- 
ments are used to compare the charmed meson and 
baryon behaviour with Sibyll 2.2c. 

The two bubble chamber experiments LEBC-EHS 
[13 at Eiab = 400 GeV and LEBC-MPS 18] at 
Eiab = 800 GeV used p-p interactions to search for 
production of all D particles. Figure[T]compares the D 
production of Sibyll 2.2c to xp and p^ distributions 
of the differential cross sections. All the comparisons 
show good agreement with data. 

The E769 experiment [19] used p, tt^, projec- 
tiles on a variety of nuclear targets to obtain the xp 
and p^ behaviour of the D mesons, where the dif- 
ferential cross sections are given per target nucleon. 
Within the perturbative QCD treatment of heavy par- 
ticle production, the incident hadron sees the nucleus 
as individual nucleons. As a consequence, the pro- 



duction cross section is expected to scale with the nu- 
clear mass A. Such behaviour has been verified by 
the E769 and WA82 ^ experiments. Nuclear scal- 
ing of the target has been upheld in the simulations, 
and a proton target has been used to compare with the 
data. Figure [5] shows the D production of Sibyll 2.2c 
plotted against the E769 data. The overall agreement 
with data is good, although the simulated p^ spectra 
are somewhat softer than data. 

The SELEX experiment [IB] used tt~ and pro- 
jectiles at Eiab = 615 GeV and p projectile at Eiab = 
540 GeV, on a variety of nuclear targets to obtain the 
distribution of charmed baryon. The leading quark 
effect in charm hadroproduction of baryons can be 
tested, as the lack of anti-quarks in the baryon beams 
and p favours quark over anti-quark and an asym- 
metry between A+ and A~ production is expected. 
SELEX observed such an asymmetry with the baryon 
beams and reported a relatively flat and hard x p dis- 
tribution, fitting to the parameterisation of (l—xp)'^'^ . 
As only the number of events of the distribution of a;^ 
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and pip are given with no absolute cross section mea- 
surement, the normalisation suggested by Ref. [2l| 
has been adopted, where the SELEX data are multi- 
plied by 7/^b/nucleon. In the simulation, 7r~ and p 
beams are used against nitrogen nuclei target. Com- 
parison of SiBYLL 2.2c with the data is shown in Fig.[3l 
For the proton beam, the shape of both x f and p^ for 
A+ distribution are in agreement though they appear 
to be over-produced using the current normalisation. 
For the ■k~ beam, the are underproduced with a 
slightly harder p^ distribution. Due to the lack of 
cross section information, we are less concerned with 
the overall scaling. The simulation agrees with the 
general trend of the data. 

IV. SUMMARY 

A new version of Sibyll 2.2c has been presented. 
This version is updated from the existing version 2.1 
to include production of charmed quarks. In addition, 
several minor changes have been made to give a bet- 
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FIG. 2: Differential cross section distributions of xp (upper panels) and (lower panels) for D production per target 
nucleon for Sibyll 2.2c. The interactions p-p (left panels), vr-p (centre panels), and K-p (right panels) are compared 
with the data from E769. 




FIG. 3: Differential cross section distribution of xp (upper panels) and (lower panels) with p (left panels) and tt 
(right panels) projectiles. The A J ( A^) production from Sibyll 2.2c shown in red solid (blue dotted) lines are compared 
with data from SELEX. 
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